An Optimal Pairs- Trading Rule* 

Qingshuo Song''' Qing Zhang* 

February 26, 2013 



Abstract 

This paper is concerned with a pairs trading rule. The idea is to monitor 
two historically correlated securities. When divergence is underway, i.e., one stock 
moves up while the other moves down, a pairs trade is entered which consists 
of a pair to short the outperforming stock and to long the underperforming one. 
Such a strategy bets the "spread" between the two would eventually converge. 
In this paper, a difference of the pair is governed by a mean-reverting model. 
The objective is to trade the pair so as to maximize an overall return. A fixed 
commission cost is charged with each transaction. In addition, a stop-loss limit 
is imposed as a state constraint. The associated HJB equations (quasi-variational 
inequalities) are used to characterize the value functions. It is shown that the 
solution to the optimal stopping problem can be obtained by solving a number of 
quasi-algebraic equations. We provide a set of sufficient conditions in terms of a 
verification theorem. Numerical examples are reported to demonstrate the results. 
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1 Introduction 



This paper is concerned with pairs trading. The idea is to identify and monitor a pair 
of historically correlated stocks. When the two stock prices diverge (one stock moves up 
while the other moves down), the pairs trade would be triggered: to short the stronger 
stock and to long the weaker one betting the eventual convergence of the prices. The 
pairs trading was first developed by Bamberger and followed by Tartaglia's quantitative 
group at Morgan Stanley in the 1980s. A major advantage of pairs trading is its 'market 
neutral' nature in the sense that it can be profitable under any market conditions. There 
are many good discussions in connection with the cause of the divergence and subse- 
quent convergence. We refer the reader to the paper by Gatev et al. [8], the book by 
Vidyamurthy [16], and references therein. 

In pairs trading, it is important to determine when to initiate a pairs trade (i.e., how 
much divergence is sufficient to trigger a trade) and when to close the position (when to 
lock in profits if the stocks perform as expected or when to cut losses if the trade goes 
sour). It is the purpose of this paper to focus on the mathematics of pairs trading. In 
particular, we consider the case when a difference of a pair satisfies a mean reversion 
model, follow a dynamic programming approach to determine these key thresholds, and 
establish their optimality. 

Mean-reversion models are often used in financial markets to capture price movements 
that have the tendency to move towards an "equilibrium" level. There are many studies 
in connection with mean reversion stock returns; see e.g., Cowles and Jones [3]) Fama and 
French [6], and Gallagher and Taylor [7] among others. In addition to stock markets, 
mean-reversion models are also used to characterize stochastic volatility (Hafner and 
Herwartz [10]) and asset prices in energy markets (see Blanco and Soronow [1]. See also 
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related results in option pricing with a mean-reversion asset by Bos, Ware and Pavlov 
[2]- 

Mathematical trading rules have been studied for many years. For example, Zhang 
[17] considered a selling rule determined by two threshold levels, a target price and a 
stop-loss limit. In [17], such optimal threshold levels are obtained by solving a set of two- 
point boundary value problems. Guo and Zhang [9] studied the optimal selling rule under 
a model with switching Geometric Brownian motion. Using a smooth-fit technique, they 
obtained the optimal threshold levels by solving a set of algebraic equations. These papers 
are concerned with the selling side of trading in which the underlying price models are of 
GBM type. Recently, Dai et al. [4] developed a trend following rule based on a conditional 
probability indicator. They showed that the optimal trading rule can be determined 
by two threshold curves which can be obtained by solving the associated Hamilton- 
Jacobi-Bellman (HJB) equations. Similar idea was developed following a confidence 
interval approach by Iwarere and Barmish [12]. In addition, Merhi and Zervos [14] studied 
an investment capacity expansion/reduction problem following a dynamic programming 
approach under a geometric Brownian motion market model. Similar problem under a 
more general market model was treated by L0kka and Zervos [13]. In connection with 
mean reversion trading, Zhang and Zhang [18] obtained a buy-low and sell-high policy 
by charactering the 'low' and 'high' levels in terms of the mean reversion parameters. 

Despite much progress in various mathematical trading rules, an important issue 
hasn't received much attention in the literature: How to cut losses and how to trade with 
cutting losses. In practice, there are many scenarios that cutting losses may arise. A 
typical one is margin call. When the pairs position is undergoing heavy losses, a margin 
call may be enforced to close part or the entire position. In addition, a pairs trader 
may determine a fixed stop-loss level from a pure money management consideration. 
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Furthermore, a historically correlated pairs may cease to be correlated at some point. 
For example, acquisition (or bankruptcy) of one stock in the pairs position. In this case, 
it is necessary to modify the trading rule to accommodate a pre-determined stop-loss 
level. From a control theoretical point of view, adding a stop-loss level is amount to 
impose a hard state constraint. This typically poses substantial difficulties in solving the 
problem. A major portion of this paper is devoted to address this important issue. 

In this paper, we consider an optimal pairs trading rule in which a pairs (long-short) 
position consists of a long position of one stock and a short position of the other. The 
state process Z t is defined as a difference of the stock prices. The objective is to initiate 
(buy) and close (sell) the pairs positions sequentially to maximize a discounted payoff 
function. A fixed (commission or slippage) cost will be imposed to each transaction. As 
in [18], we study the problem following a dynamic programming approach and establish 
the associated HJB equations for the value functions. We show that the corresponding 
optimal stopping times can be determined by three threshold levels xq, x±, and x 2 . These 
key levels can be obtained by solving a set of algebraic like equations. We show that 
the optimal pairs trading rule can be given in terms of two intervals: l\ = [x ,Xi] and 
J 2 = (M, x 2 ). Here M is the given stop-loss level (e.g., as the consequence of a margin 
call) and I\ is contained in J 2 . The idea to initiate a trade whenever Z t enters I\ and 
hold the position till Z t exits J 2 . In addition, we provide a set of sufficient conditions 
that guarantee the optimality of our pairs trading rule. We also examine the dependence 
of these threshold levels on various parameters in a numerical example. Finally, we 
demonstrate how to implement the results using a pair of stocks and their historical 
prices. 

This paper is organized as follows. In §2, we formulate the pairs trading problem 
under consideration. In §3, we study properties of the value functions, the associate 
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HJB equations, and their solutions. In §4, we provide a set of sufficient conditions that 
guarantee the optimality of our trading rule. A numerical example is given in §5. The 
paper is concluded in §6. 

2 Problem Formulation 

Let X] and X 2 denote the prices of a pair of correlated stocks X 1 and X 2 , respectively. 
The corresponding pairs position consists of a long position in stock X 1 and short position 
in stock X 2 . For simplicity, we include one share of X 1 and K shares of X 2 in the pairs 
position. Here K Q is a given positive number. The price of the position is given by 
Z t = X} — KqX 2 . We assume that Z t is a mean-reverting (Ornstein-Uhlenbeck) process 
governed by 

dZ t = a{b - Z t )dt + adW t , Z = x, (1) 

where a > is the rate of reversion, b the equilibrium level, a > the volatility, and W t 
a standard Brownian motion. 

In this paper, the notation X 1 , i = 1,2, are reserved for the underlying stocks and 
Z the corresponding pairs position. One share long in Z means the combination of one 
share long position in X 1 and K shares of short position in X 2 . Similarly, for % — 1, 2,, 
X\ represents the price of stock X % and Z t the value of the pairs position at time t. Note 
that Z t is allowed to be negative in this paper. 

In addition, we impose a state constraint and require Z t > M. Here M is a given 
constant and it represents a stop-loss level. It is common in practice to limit losses to 
an acceptable level to account for unforeseeable events in the marketplace. A stop-loss 
limit is often enforced as part of money management. It can also be associated with a 
margin call due to substantial losses. 



5 



To accommodate such state constraint in our model, let tm denote the exit time of 
Z t from (M, oo), i.e., r M = inf{t : Z t g (M, oo)}. 
Let 

< r\ < rf < t\ < r| < ■ ■ ■ < t m (2) 

denote a sequence of stopping times. A buying decision is made at and a selling 
decision at r*, n — 1, 2, 

We consider the case that the net position at any time can be either long (with one 
share of Z) or flat (no stock position of either X 1 or X 2 ). Let i = 0, 1 denote the initial 
net position. If initially the net position is long (i = 1), then one should sell Z before 
acquiring any future shares. The corresponding sequence of stopping times is denoted 
by A x = (t* , r 2 > r 3> • • •)• Likewise, if initially the net position is flat {i = 0), then 
one should start to buy a share of Z. The corresponding sequence of stopping times is 
denoted by A = (rf , rf, r|, r|, . . .). 

Let K > denote the fixed transaction cost (e.g., slippage and/or commission) asso- 
ciated with buying or selling of Z. Given the initial state Z = x and initial net position 
i — 0,1, and the decision sequences, A and A 1; the corresponding reward functions 



Ji(x,hi) 



El £ [e-^(Z T s -K)- e-ri{Z Ti + K)] I {<<tm} , if i = 0, 



E\e-^{Z T s-K) (3) 



+ ^[e- pT HZ^-K)-e-^(Z T , + K) 



n=2 

where p > is a given discount factor. 

In this paper, given random variables £ n , the term E Z)J£=i is interpreted as 

N 

limsup£^£„. 
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In the reward function J { , a buying decision has to be made before Z t reaches M. 
When t = t m (or Z t = M), only a selling can be done if % — 1. 

For i = 0,1, let Vi(x) denote the value functions with the initial state Z = x and 
initial net positions i — 0,1. That is, 

Vi(x) = supJi(x,Ai). (4) 

A; 

Note that 

V Q (M) — and Vi(M) — M — K. (5) 
These give the boundary conditions. 

Remark 1. Note that we allow the equalities in (2), i.e., one can buy and sell simul- 
taneously. Nevertheless, owing to the existence of positive transactions cost K, any 
simultaneous buying and selling are automatically ruled out by our optimality condi- 
tions. 

We also imposed the conditions < t m and < t m , n = 1,2, If one has a 

share position of Z and = tm for some n, then one has to sell the share to cut losses. 
On the other hand, if r h n = t m , then one should not buy because she has to sell it right 
away, which only cause the round trip transaction fees. 

Remark 2. Recall that in this paper the stock (pair) price is given by Z t . In [18], a 
percentage slippage cost is required and the stock price is given by S t = e Zt . Suppose 
K percentage is added to a buying order. Then the total cost is given by S t (l + K) = 
e Zt (l + K). Its natural logarithm equals approximately Z t + K, which matches the cost 
structure in this paper. 

Remark 3. In addition, we only consider the 'long' side trading in this paper. Actually, 
one can trade by simply reversing the trading rule obtained in this paper. For example, 
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if b = 0, then we can trade both Z t and (—Z t ) simultaneously because they satisfy the 
same system equation (1). 

Remark 4. The optimal stopping problem considered in this paper can be generalized 
to treat similar problems in related fields (e.g., the energy market). We refer the reader 
to Hamadene and Zhang [11] and references therein for additional applications. 

Example 1. Typically a highly correlated pair can be found from the same industry 
sector. In this example, we choose Wal-Mart Stores Inc. (WMT) and Target Corp. 
(TGT). Both companies are from the retail industry and they have shared similar dips 
and highs. If the price of WMT were to go up a large amount while TGT stayed the same, 
a pairs trader would buy TGT and sell short WMT betting on the convergence of their 
prices. In Figure 1, the 'normalized' (dividing each price by its long term moving average) 
difference of WMT and TGT is plotted. In addition, the data (1992-2012) is divided into 
two sections. The first section (1992-2000) is used to calibrate the model and the second 
section (2001-2012) to backtest the performance of our results. Our construction of Z t 
determines that the equilibrium level 6 = 0. By measuring the standard derivation of 
Z t , we obtain the historical volatility a = 0.56. Finally, following the traditional least 
squares method, we obtain a = 1.00. 

3 Properties of the Value Functions 

In this section, we establish various bounds for the value functions and solve the associ- 
ated HJB equations. 

First, note that the sequence A = (rf, r-f, t%, t|, . . .) can be regarded as a combination 
of a buy at t\ and then followed by the sequence of stopping times A 1 = (r-f, r|, r|, r|, . . .). 
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1 - 1992/05/08 — 2000/12/29 



2001/01/02 — 2012/12/31 



Figure 1: WMT and TGT (1992-2012) 

In view of this, we have, for x > M, 
V (x) > J (x,A ) 

{QO 
e-rt(Z T s - K)I {Tl<TM} + £ [e-»<{Z T s - K) - e"^(^ + K)] 7 {t , <Tm} 
n=2 

-Ee-ri(Z Tbi +K)I {Tbi<TM} 
= J,(x, A,) - Ee-^{Z rl + K)I {rl<TM y 
In particular, setting r\ = and taking supremum over A l5 we obtain the inequality 



V (x) > V^-x-K. 



(6) 



Similarly, we can show, for x > M, that 



Vi{x) > V (x)+x-K. 



(7) 



Clearly, in view of the boundary conditions (5) these two inequalities hold for x = M. 
Next, we establish lower and upper bounds for Vi(x). 
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Lemma 1. The following inequalities hold: 

< V (x) < C , 

x-K < V^x) <x + K + C , 

for all x G [M, oo), where C = (p + a)\M\/p. 

Proof. Note that the lower bounds for Vi(x), {i = 0, 1), follow from their definitions. 
In addition, if C is an upper bound for Vo(x), then the upper bound for Vi(x) follows 
from the inequality in (6). It remains to show the upper bound for Vq. Recall that 
T~n < < tm- Therefore, we have 

E (W Ti - W*) I {T , <TM} = E (W T s - W T >) - E (W T s - W T >) I {Tbn=TM} = 0. 

Recall also that Z t > M for all t < t m . Using Dynkin's formula, we have, for each 
n= 1,2,..., 

£(e^^-e-^)/ tt<TM} 

= E (/f e ' pt (~(P + a ) Z t) dt ) kri<r M } + E (a(W T s - W T >)) J K<TM} 
<(p + a)\M\E^e-^dtjl {Th<TM} 
< (p + a)\M\E / e~ pt dt. 

It follows from the definition of Jo(x, A ) that 

oo 

J (x, Ao) < Y.E {e~^Z T s - e-riz T >) 7 {rf<TM} 

n=l 

OO r. T h 

< (p + a)\M\Y,E n e- pt dt 

n=l Jr n 

PCX) 

< (p + a)\M\ e~ pt dt 

(p + a)\M\ °_ 
— — <^o- 



(8) 
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pvo — Avo = vq = v\ — x—K pvQ — Avo = 

l 1 1 - 

M xq X\ 

l 1 - 

pv\ — Av\ = I) v\ = vq+x — K 

Figure 2: Continuation Regions 
This implies V (x) < C . □ 

Let A denote the generator of Z t , i.e., 

A = a(o — x)— — h 



(9) 



dx 2 dx 2 

Formally, the associated HJB equations should have the form: 

min ^pvo(x) — Avo(x), vq(x) — v\{x) + x + = 0, 
min j/w i(x) — Av±(x), Vi(x) — vq(x) — x + = 0, 

for x G (M, oo), with the boundary conditions v (M) = and V\{M) = M — K. 

If i — 0, then one should only buy when the price is low (say less than or equal to 
xi). In this case, Vq{x) — V\{x) — x — K. The corresponding continuation region (given 
by pvo(x) — .4^0(2;) = 0) should include (xi,oo). In addition, one should not establish 
any new position if Z t is close to the stop-loss level M. In view of this, the continuation 
region should also include (M,xq) for some xq < x\. On the other hand, if % — 1, then 
one should only sell when the price is high (greater than or equal to x-i > xi), which 
implies Vi(x) = v (x) +x — K and the continuation region (given by pv\{x) —Avi(x) = 0) 
should be (M, x?). These continuation regions are highlighted in Figure 2. 

To solve the HJB equations in (9), we first solve the equations pvi(x) — Avi(x) = 
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with i = 0, 1 on their continuation regions. Let 



faix) = / V (t)e-< b -^dt, 
Jo 

^(x) = / v {t)e K(b - x)t dt, 
Jo 



where r)(t) = t^^^ 1 exp (— 1 2 /2) and k = y/2a/a. Then the general solution (see Eloe et 
al. [5]) is given by A\$\{x) + A 2 (j) 2 (x), f° r some constants A\ and A®. 

First, consider the interval (x±, oo) and suppose the solution is given by + 
A 2 (p 2 (x), for some A\ and A 2 . Recall the upper bound for Vo(x) in Lemma 1, t>o(oo) 
should be bounded above. This implies that, A\ — and vq{x) = A 2 <p 2 (x) on (xi,oo). 
Let B-l, B 2 , Ci, and C 2 be constants such that v (x) = Bi<j>i(x) + B 2 (f) 2 (x) on (M, x ) 
and = Ci0i(x) + C 2 02(^) on (x 2 , oo). 

It is easy to see that these functions are twice continuously differentiable on their 
continuation regions. We follow the smooth-fit method which requires the solutions to be 
continuously differentiable. In particular, it requires v to be continuously differentiable 
at x . Therefore, 



Bi^i(x ) + B 2 (p 2 (x ) = Ci0i(x o ) + C 2 (/) 2 (xo) -x - K, 



(10) 



[ BMxo) + B 2 <j>' 2 (x ) = CMxo) + C 2 <f>' 2 (x ) - 1. 
Similarly, the smooth-fit conditions at x\ and x 2 yield 



^202(Zl) = C , i0i(xi) + C 2 2 (^i) - Xi 



K, 



(11) 



< 



A 2 <i>' 2 (x 1 ) = CMxJ + C 2 2 (^i) - 1, 



and 



Ci(f> 1 (x 2 ) + C 2 (fi 2 (x 2 ) = A 2 (p2(x2) + x 2 



K, 



(12) 



_ Ci0i(^2) + C 2 <t>' 2 {x 2 ) = A 2 <t>' 2 {x 2 ) + 1. 



Finally, the boundary conditions at x = M lead to 




(13) 
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For simplicity in notation, let 



f M x ) </> 2 (x) ^ 



V <P'l( X ) 02 ( X ) J 

Note that the determinant of <&(x) is given by 



-K 



r v (t)e- K{b - x)t dt rtr]{t)e K{b - x)t dt+ f°° t^e^-^dt f° V (t)e K{b 
o Jo Jo Jo 



which is less than zero for all x. Therefore, <&(x) is invertible for all x. 
Also, let 



R( x ) = 



, P 1 (x) = <S>-\x) 



4+^ 



v 1 J 



( 



, P 2 (x) = 



x- 



Rewrite the equations (10)-(13) in terms of these vectors. We have 



\B 2 ) 



A 2 R{ Xl ) 



\ ° 2 ) 



-PiM, 



V ° 2 ) 



-Pi(xi), 



V c * J 

A 2 R{x 2 )+P 2 {x 2 ), 



and 



(MM), MM)) 



(MM), MM)) 



f B1 \ 

K B 2J 
( Cl \ 

\ C * ) 



o, 



= M-K. 
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Multiplying both sides of (14) from the left by (0i(M), <f> 2 (M)) and using (17), we 
have 

(MM),MM))Pi(x ) = M-K. (18) 
Combining (15) and (16) and eliminating (Ci,C 2 ) T , we obtain 

A 2 (R( Xl ) - R(x 2 )) = P 2 (x 2 ) - P l (x l ). (19) 

Also, multiplying both sides of (16) from the left by (0i(M), (f) 2 (M)) yields 

M — K — A 2 ((f>i(M), 4> 2 {M))R{x 2 ) + (MM), 4> 2 {M))P 2 {x 2 ). (20) 

It is easy to check that 

(MM),MM))R(x 2 ) = 2 (M) det$(:r 2 ) ^ 0. 

This leads to 

A M-K-(MM),MM))P 2 (x 2 ) 

(MM),MM))R(x 2 ) ■ { ] 

Finally, substitute this into (19) to obtain 

/ / n / (M-K- (6 1 (M),(j) 2 (M))P 2 (x 2 )\ 

Solving equations (18) and (22), we can obtain the triple (xo, xi, x 2 ). Then solving the 
equations (14), (15), and (21), to obtain A 2 , (B 1: B 2 ), and (C 1: C 2 ). 

We need additional conditions for x x and x 2 . Note that Vi(x) has to satisfy the 
following inequalities for being solutions to the HJB equations (9): 

pv (x) - Av (x) > 0, 
pvi(x) — Avi(x) > 0, 

(23) 

v (x) > Vi(x) — x — K, 
Vi(x) > Vq(x) + x — K, 
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for all x > M. Next, we examine each of these inequalities on intervals (M,x ), (x ,Xi), 
(x 1 ,x 2 ), and (x 2 , oo). 

First, on (M, x ), the top two inequalities in (23) become equalities. We only need 
the last two inequalities to hold. Therefore, we have 

x - K < vi(x) - v (x) < x + K on {M,x ). (24) 

Then, on (xo,xi), note that vo(x) = v\{x) — x — K implies vi(x) > vo(x) +x — K. We 
only need pvo(x) — Av (x) > 0. Again, using vo(x) = vi(x)—x — K and pv i(x) — Av\(x) = 
on this interval, we have 

pv (x) — Av (x) = p(vi(x) — x — K) — A(vi(x) — x — K) 
= p(—x — K) — A(—x — K) 
= — (p + a)x — pK + ab. 
In view of this, pvo(x) — Avq{x) > on (xq,x\) is equivalent to 

Xl < ah - pK . (25) 
~ p + a 

Similarly, on (xi,x 2 ), we only need the inequalities 

x - K < vi(x) - v (x) <x + K. (26) 
Finally, on (x 2 , oo), we only require 

ab + pK . 

x 2 > — . 27 

p + a 

Note that the inequalities in (24) and (26) are equivalent to the following inequalities, 

\(C 1 -B 1 )(f> 1 (x) + (C 2 -B 2 )<h(x)-x\<K on(M,x ), 

(28) 

\d<f>i(x) + (C 2 - A 2 )<f) 2 (x) - x\ < K on (xi,x 2 ), 

respectively. 

In what follows, we show that the triple (x ,Xi,x 2 ) satisfying these conditions leads 
to the optimal stopping rules. 
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4 A Verification Theorem 



V (x) 



In this section, we give a verification theorem to show that the solution Vi(x), % — 0, 1, of 
equation (9) are equal to the value functions Vi(x), i = 0,1, respectively, and sequences 
of optimal stopping times can be constructed from the triple (x , x±, x 2 ). 

Theorem 1. Let (x ,xi,x 2 ) be a solution to (18) and (22) and satisfy 

ah — pK ab + pK 

X\ < and x 2 > . 

p + a p+ a 

Let A 2 , B\, B 2 , C\, and C 2 be constants given by (14), (16), and (21) satisfying the 
inequalities in (28). 
Let 

B^x) + B 2( f> 2 (x) ifxe[M,x ), 
Ci<pi{x) + C 2 (f) 2 (x) - x - K if x G [x ,xi), 
A 2 <p 2 (x) ifxe [xi, oo), 

CMx) + C 2 <p 2 (x) ifxe [M,x 2 ), 
A 2 (f) 2 (x) + x - K if x G [x 2 , oo). 
Assume v (x) > 0. Then, Vi(x) = Vi(x), i — 0,1. Moreover, if initially i = 0, let 

a* / b* s* b* s* \ 

Y H) ~ \ T 1 ) r l ? r 2 ' r 2 J • • •)■> 

such that the stopping times r b * = inf{t > : Z t G [x ,Xi]} A t m , t£* = inf{£ > r b * : 
Z t (M, x 2 )} A r M; and r b * +1 = inf{t > r** : Z t G [x , Xi]} A r M for n > 1. Similarly, if 
initially i — 1, let 

such that t(* = inf{t > : Z t £ (M, x 2 )} A t m , r b * = M{t > r** x : G [x , Xi]} A r M , 
and r^* = inf{t > r b * : Z t G" (M, x 2 )} A t m for n>2. Then Aq and A* are optimal. 



Vi(x) 
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Proof. We divide the proof into two steps. In the first step, we show that Vi(x) > Ji(x, Aj) 
for all Aj. Then in the second step, we prove that Vi(x) = Jj(x, A*), which implies 
Vi(x) = Vi(x) and A* is optimal. 

Let Io = (M,xo) U (xo,xi) U (xi,oo) and l\ = (M, x 2 ) U (x2,oo). It is easy to see 
that v G C 2 (I ), vi E C 2 (Ii), and both v and v\ are in C 1 ([M, oo)). In addition, 
they satisfy the quasi-variational inequalities in (9), i.e., pvi{x) — Avi(x) > 0, i — 0, 1, 
whenever they are twice continuously differentiable. Using these inequalities, Dynkin's 
formula, and Fatou's lemma as in 0ksendal [15, p. 226], we have, for any stopping times 
< 9i < 02 < tm, a.s., 



Ee-&Vi(Z ei ) > Ee'^v t (X e2 ), 
Ee-^ Vi (Z ei )I { o 1<TM} > Ee-^ Vi (X e2 )I {ei<TMh 

for i — 0,1. Given A = (r-f , r-f , r\, r|, . . .), using (6) and v (M) = 0, we have 

v (x)> Ee- pT 'v (Z T >) 

= Ee-ri Vo (Z T >)I {Tbi<TM} 

> Ee-ri ( Vl (Z Tbi ) - {Z rl + K)) / {rf<TM} 

= Ee-ri Vl (Z T >)I {Tbi<TM} - Ee-ri(Z Tbi + K)I {rb<TM} . 

It follows again from (29) and then (7) that 

v (x) > Ee-rt Vl (Z Tl )I {Tbi<TM} - Ee~ri{Z Tb + K)I {rb<TM} 

> Ee-^ {vo{ Z T s) + Z T s - K)I {rb<TM} - Ee-ri{Z rb + K)I {rb<TM} 

= Ee-^v (Z T s)I {Tb<TM} + E[e-^(Z T s — K) — e~^{Z rb + K)]l {rb<TM} 

= Ee-rtv (Z T s) +E[e-rt(Z T s — K) — e~^{Z rb + K)\l {rb<TM y 

Note that 

Ee-^v (Z Tl ) > Ee-^v (Z Tb ) = Ee-^v {Z Tb )I {rb<TM} . 
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(29) 



Similarly, we have 



Ee- pT ?v (Z T s) > Ee- pT iv (Z Ti ) + E e~ p ^(Z T s - K) - e~ pT *{Z T b + K) 



Repeat this process and note that v (x) > to obtain 



h-4<TuY (3°) 



N 

v (x) >EY: [e- pT HZr* -K)- e-ri(Zr> + K) 



71=1 



l {rk<T M }- 



Sending N — > oo to obtain vo(x) > Jq(x,A ) for all A . Therefore, vo(x) > Vo(x). 
Similarly, using (30), we can show that 

vi(x) > Ee- pT ? Vl (Z T s) 

> Ee- pT ? (v (Z T s) + Z T s - K) 

= Ee- pT f (Z T s -K) + Ee- pT ?v (Z Tl ) 

> ■■■ 



N 

= Ee~^ (Z T s -K)+EY: [e-« (Z T * - K) - e~ p < (Z T » + K) 



n=2 



1 {t*<t m }- 



It follows that Vi(x) > Vi(x). 

Next, we establish the equalities. Define t\* = inf{i > : Z t e [x ,Xi]} A t m . Note 
that r M < oo, a.s. (see [18, Lemma 6]). Therefore, rf* < oo, a.s. Using again Dynkin's 
formula, we have 

v (x) = Ee- pT i*v (Z T b*) 

= Ee-ri*v (Z Tt )I {Tbi<TM} 

= Ee-rf* (y^) - (Z Tt + K)) 7 {rf<TM} 

= Ee-ri* Vl (Z Tt )I {Tbi<TM} - Ee-rt\Z Tt + K)I {t1<tm} . 
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Let rf* = mi{t > t\* : X t = x 2 } A r M . Then, rf* < oo, a.s. We have also 



Ee-^ Vl (Z Tbi *)I {Tbi<TM} 

= Ee^r Vl{Zrr )I {Tbi<TM} 

= Ee'rt* (v (Z rr ) + (Z rr - K)) I {Tb<TM} 

= Ee-^t Vo{ Z Tr )I {Tbi<TM} + Ee-^\ Zrr - K)I {rb<TM} 

= Ee-^*v Q {Z TV ) + Ee-^\z TV - K)I {rb<TM} . 

It follows that 



v (x) = Ee- pT "v (Z rr ) + E 



K) 



e-^{Z Tb , 



K) 



! {^<r M }- 



Continue this way to obtain 

v (x) = Ee-^MZrj.) + £ £ [e- pT "(Z rr - K) - e~^{Z^ + K)]l {rt :] . 



N 



n=l 



Similarly, we can show 

Vl (x)= Ee-rt' Vl (Z rr ) 

= Ee-^(v (Z rr ) + Z rr -K) 

= Ee~ pT "v (Z Tr ) + Ee~ pT °* (Z T .. - K) 

= Ee~ pT N vo(Z T s* ) + Ee- pT ?* (Z T s* - K) 

N 

+E ]T [e~ pT " {Z r s, -K)- e~ p <* (Z rb * + K) 

n=2 

Recall that P(r M < oo) = 1. This implies lini7v->oo r^* = r M , a.s. Recall also that 
v (M) = 0. It follows that Ee~ pT " 'v (Z T **) -»■ 0. This completes the proof. □ 

5 A Numerical Example 

In this section, we use the parameters of the WMT-TGT example, i.e., 
a = 1.0, 6 = 0, a = 0.56, p = 0.10, K = 0.001. 
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Solving the equations (18) and (22) gives the triple (x , Xi, x 2 ) = (—0.142, —0.077, 0.077). 
Next, we vary one of the parameters at a time and examine the dependence of the triple 
(x ,Xi,x 2 ) on these parameters. 

Dependence of (^0,^1,^2) on parameters 

First we consider the triple (xo,xi,x 2 ) associated with varying a. A larger a implies 
larger pulling rate back to the equilibrium 6 = 0. It can be seen in Table 1 that the lower 
buying level xq decreases as a gets bigger. Also the higher buying level x\ increases in a. 
These lead to larger buying interval [xoj^i] resulting greater buying opportunities. The 
selling level x 2 decreases which suggests one should take profit sooner as a gets bigger 
because the potential of going higher becomes smaller. In addition, the interval (xi,x 2 ) 
is symmetric about b = 0. 



a 


0.60 


0.80 


1.00 


1.20 


1.40 


x 


-0.124 


-0.135 


-0.142 


-0.147 


-0.151 


Xl 


-0.089 


-0.083 


-0.077 


-0.073 


-0.069 


x 2 


0.089 


0.083 


0.077 


0.073 


0.069 



Table 1. (xo,x±,x 2 ) with varying a. 



In Table 2, we vary the volatility a. The volatility is the source forcing the price to go 
away from its equilibrium. The large the a, the further the price fluctuates. As a result, 
every element in the triple (x ,Xi,x 2 ) moves along the opposite direction as a increases 
resulting a smaller buying interval [xq,xi] and a higher profit target x 2 . 
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a 


0.36 


0.46 


0.56 


0.66 


0.76 


Xo 


-0.164 


-0.153 


-0.142 


-0.130 


-0.117 


Xl 


-0.057 


-0.067 


-0.077 


-0.086 


-0.095 


X2 


0.057 


0.067 


0.077 


0.086 


0.095 



Table 2. (xq,xi,X2) with varying a. 



Next, we vary the discount rate p. Larger p means quicker profits. This is confirmed 
in Table 3. It shows that larger p leads to a smaller x , a slightly larger x±, and a slightly 
smaller x-i- This means more buying opportunities and quicker profit taking. 



p 


0.06 


0.08 


0.10 


0.12 


0.14 


x 


-0.1412 


-0.1416 


-0.1420 


-0.1426 


-0.1430 




-0.078 


-0.078 


-0.077 


-0.077 


-0.076 


X2 


0.078 


0.078 


0.077 


0.077 


0.076 



Table 3. (xq,x±,X2) with varying p. 



Finally, we examine the dependence on the stop-loss level M. Clearly, a smaller M 
is associated with a larger loss when it goes wrong. In Table 4, the lower buying level xo 
decreases in M. On the other hand, the buying-selling interval (x\, X2) is not as sensitive 
to variations in M. 



M 


-0.16 


-0.18 


-0.20 


-0.22 


-0.24 


xo 


-0.091 


-0.118 


-0.142 


-0.166 


-0.189 


Xl 


-0.077 


-0.078 


-0.077 


-0.077 


-0.077 


X2 


0.077 


0.078 


0.077 


0.077 


0.077 



Table 4. (xq, x\, X2) with varying M. 
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Backtesting (WMT-TGT) 



We backtest the pairs trading rule using the stock prices of WMT and TGT from 2001 
to 2012. Let Xl be the WMT stock divide by its 1000 day moving average and 
the TGT stock by its same period moving average. We take Z t = X\ — Xf. Using the 
parameters obtained in Example 1 based on the historical prices from 1992 to 2000, we 
found the triple (x ,xi,x 2 ) = (—0.142, —0.077,0.077). A pairs trading is triggered when 
Z t gets inside the buying interval [x , Xi\. The position is closed when Z t exits the interval 
(M, x 2 ). Initially, we allocate trading the capital $100K. When the first long signal is 
triggered, buy $50K WMT stocks and short the same amount TGT. Close the position 
either when Z t reaches the target x 2 or when it drops below the stop-loss level M. Such 
half-and-half capital allocation between long and short applies to all trades. In addition, 
each pairs transaction is charged $5 commission fee. Furthermore, two variations from 
the assumptions prescribed in Theorem 1 in our 'actual' trading: (a) After the stop-loss 
level M is reached, the trading continues and a buying order is entered when Z t goes 
back to the trading range; (b) All available capital will be used (half long and half short) 
for trading rather than following the 'single' share rule, 

In Figure 3, the corresponding Z t , the threshold triple, and the corresponding equity 
curve are plotted. There are total 8 trades and the end balance is $126.602K. 

Note that Z t is symmetric, i.e., (—Z t ) satisfies the same equation (1). Naturally, one 
can reverse the pair and trade (—Z t ) the same way. The reversed Z t and equity curve is 
given in Figure 4. Such trade leads to the end balance $114.935K. Note that both types 
of trades have no overlap, i.e., they do not compete for the same capital. The grand total 
profit is $41547 which is a %41.54 gain. 
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Equity curve (X 100K) : 2001/01/02 — 2012/12/31 
(Long WMT and short TGT ) 



I AmA f '\ ft, w i 
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Figure 3: Threshold levels and the equity curve 



Equity curve (X 100K) : 2001/01/02 — 2012/12/31 
(Long TGT and short WMT) 
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Figure 4: Threshold levels and the equity curve 
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The main advantage of pairs trading is its risk neutral nature, i.e., it can be profitable 
regardless the general market condition. In addition, there are only 2x8 trades leaving 
the capital in cash most of the time. This is desirable because the cash sitting in the 
account can be used for other types of shorter term trading in between, at least drawing 
interest over time. 

Finally, the choice of stop-loss level M can depend on many factors including the 
trader's risk tolerance level and margin requirements. Our choice M = —0.2 corresponds 
to a %10 loss when WMT drops %10 and TGT stays the same. 

6 Conclusion 

In this paper, we have studied the pairs trading problem following a mean reversion ap- 
proach and obtained a closed-form solution under reasonable conditions. Much attention 
was given to the trading rule with loss cutting, which is an important component of 
money management. 

A simple real market (WMT-TGT) example was considered. It would be interesting 
to examine how the method works for a larger selection of pairs of correlated stocks. 
Some practical considerations can be found in the book by Vidyamurthy [16]. 
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